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TOEPLITZ AND TOEPLITZ-BLOCK-TOEPLITZ MATRICES 
AND THEIR CORRELATION WITH SYZYGIES OF POLYNOMIALS 

HOUSSAM KHALIL*, BERNARD MOURRAINt, AND MICHELLE SCHATZMAN* 

Abstract. In this paper, we re-investigate the resolution of Toeplitz systems Tu = g, from a new point of view, by 
correlating the solution of such problems with syzygies of polynomials or moving lines. We show an explicit connection 
between the generators of a Toeplitz matrix and the generators of the corresponding module of syzygies. We show that 
this module is generated by two elements of degree n and the solution of Tu = g can be reinterpreted as the remainder 
of an explicit vector depending on g, by these two generators. 

This approach extends naturally to multivariate problems and we describe for Toeplitz-block- Toeplitz matrices, the 
structure of the corresponding generators. 
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1. Introduction. Structured matrices appear in various domains, such as scientific computing, 
signal processing, . . . They usuaUy express, in a hncarize way, a problem which depends on less pa- 
rameters than the number of entries of the corresponding matrix. An important area of research 
, is devoted to the development of methods for the treatment of such matrices, which depend on the 

actual parameters involved in these matrices. 

Among well-known structured matrices, Toeplitz and Hankel structures have been intensively 
1^ ' studied [3 [6] . Nearly optimal algorithms are known for the multiplication or the resolution of linear 

, systems, for such structure. Namely, if A is a Toeplitz matrix of size n, multiplying it by a vector or 

solving a linear system with A requires 0{n) arithmetic operations (where 0{n) = O{nlog'^{n)) for 
some c > 0) |H|T2]. Such algorithms are called super-fast, in opposition with fast algorithms requiring 
0{n?) arithmetic operations. 
^ . The fundamental ingredients in these algorithms are the so-called generators [B], encoding the 

' minimal information stored in these matrices, and on which the matrix transformations are translated, 

■"sj" . The correlation with other types of structured matrices has also been well developed in the literature 

I [ini [9] , allowing to treat so efficiently other structures such as Vandermonde or Cauchy-like structures. 

Such problems are strongly connected to polynomial problems [H [1] . For instance, the product 
I of a Toeplitz matrix by a vector can be deduced from the product of two univariate polynomials, and 

■ thus can be computed efficiently by evaluation-interpolation techniques, based on FFT. The inverse 

I of a Hankel or Toeplitz matrix is connected to the Bezoutian of the polynomials associated to their 

generators. 

However, most of these methods involve univariate polynomials. So far, few investigations have 
' been pursued for the treatment of multilevel structured matrices jllj . related to multivariate prob- 

^ I lems. Such linear systems appear for instance in resultant or in residue constructions, in normal form 

computations, or more generally in multivariate polynomial algebra. We refer to [8] for a general 
description of such correlations between multi-structured matrices and multivariate polynomials. Sur- 
prisingly, they also appear in numerical scheme and preconditionners. A main challenge here is to 
devise super- fast algorithms of complexity 0{n) for the resolution of multi-structured systems of size 
n. 

In this paper, we consider block- Toeplitz matrices, where each block is a Toeplitz matrix. Such 
a structure, which is the first step to multi-level structures, is involved in many bivariate problems, 
or in numerical linear problems. We re-investigate first the resolution of Toeplitz systems T u = g, 
from a new point of view, by correlating the solution of such problems with syzygies of polynomials 
or moving lines. We show an explicit connection between the generators of a Toeplitz matrix and the 
generators of the corresponding module of syzygies. We show that this module is generated by two 
elements of degree n and the solution oiT u — g can be reinterpreted as the remainder of an explicit 
vector depending on g, by these two generators. 
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This approach extends naturaUy to multivariate problems and we describe for Toeplitz-block- 
Toeplitz matrices, the structure of the corresponding generators. In particular, we show the known 
result that the module of syzygies of k non-zero bivariate polynomials is free of rank fc — 1 , by a new 
elementary proof. 

Exploiting the properties of moving lines associated to Toeplitz matrices, we give a new point of 
view to resolve a Toeplitz-block- Toeplitz system. 

In the next section we studie the scalar Toeplitz case. In the chapter 3 we consider the Toeplitz- 
block- Toeplitz case. 

Let R = K[x]. For n S N, we denote by K[x]„ the vector space of polynomials of degree < n. 
Let L = 'K[x,x~^] be the set of Laurent polynomials in the variable x. For any polynomial p = 
J27=~mPi^^ £ L, we denote by p+ the sum of terms with positive exponents: p+ = J2^=oPi^^ 
by , the sum of terms with strictly negative exponents: p^ = ^i2-mPi have p = p'^ + P^ ■ 

For n E N, we denote by ii.„ — {lu; uj^ ~ 1} the set of roots of unity of order n. 

2. Univariate case. We begin by the univariate case and the following problem: 
PROBLEM 2.1. Given a Toeplitz matrix T = {ti-j)"'^^ £ K"^" (T = {T^jY-'Iq with T,j = t,_j) 
of size n and g = (50, • ■ ■ , 9n~i) G IK", find u = (uq, . . . , m„-i) G K" such that 

Tu^g. (2.1) 



Let E ^ {1, . . . .x"^ and IIb be the projection of R on the vector space generated by E, along 
(a;",a;"+i,...}. 

Definition 2.2. We define the following polynomials: 



n-l 

• J (X) = tiX', 

i——n+l 

2n-l . 

• f{x) = Ux' With U = \ ' ifKn 

i=0 ^ 
n-l 



ti-2n ifi>n ' 

1=0 1=0 

Notice that f = T+ + a;^ " and T{w) = f{w) if w £ 5l2„- We also have (see [5]) 

T u — g <^ IIe{T{x)u{x)) — g{x). 

For any polynomial u G K[x] of degree d, we denote it as u{x) = u{x) + x"'u{x) with deg(w) < n — 1 
and deg(M) <d — nifd>n and u = otherwise. Then, we have 

T{x) u{x) = T{x)u{x) + T{x)x"u{x) 

= IIe{T{x)u(x)) + IIe{T{x)x'^u{x)) 

+ {a^n+ix~"+^ + h a^ix-^) 

+ (a„a;" + • • • + a„+™a;"+") 
= IlEiT{x)u{x)) + IlEiT{x)x''u{x)) 

+x-''+'^A{x) + x''B{x), (2.2) 

A{x) = a_„+i -I- • • ■ -I- a-ix"-^, 

B{x)=an + --- + an+mx"'. (2.3) 



with m — max(n — 2, d — 1) 



See [8] for more details, on the correlation between structured matrices and (multivariate) polynomials. 

2.1. Moving lines and Toeplitz matrices. We consider here another problem, related to 
interesting questions in Effective Algebraic Geometry. 

PROBLEM 2.3. Given three polynomials a,b,c £ R respectively of degree < l^ < m, < n, find three 
polynomials p,q,r £ R of degree < v ~ l^ < v ~ < v ~ n, such that 

a{x)p{x) + b{x) q{x) + c{x) r{x) = 0. (2.4) 
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We denote by C{a, b, c) the set of {p, q, r) G K[x]'^ which are solutions of (|2.4p . It is a ]K[x]-niodule of 
K[a:]^. The solutions of the problem ()2.3p are C{a,b,c) nK[a;]iy_;_i x K[x]i/_,„_i x ]K[x]i/_„_i. 

Given a new polynomial d{x) G K.[x], we denote by C{a,b,c;d) the set of {p,q,r) G K[a;]^ such 

that 

a{x) p{x) + b{x) q{x) + c{x) r{x) — d{x). 



Theorem 2.4. For any non-zero vector of polynomials (a, 6, c) G K[a;]'^, the K.[x]-module C{a, b, c) 
is free of rank 2 . 

Proof. By the Hilbert's theorem, the ideal / generated by (a, 6, c) has a free resolution of length 
at most 1, that is of the form: 

^ ]K[a;]P ^ K[xf K[x\ K[x\II 0. 

As / ^ 0, for dimensional reasons, we must have p = 2. U 

Definition 2.5. A ji-base of C{a,b,c) is a basis {p,q,r), {p',q',r') of C{a,b,c), with {p,q,r) of 
minimal degree /i. 

Notice if /ii is the smallest degree of a generator and /i2 the degree of the second generator 
{p' , q' , r'), we have d = max(deg(a), deg(6), deg(c)) = /ii + /i2- Indeed, we have 

^ K[x]^^d-^,^ © K[.x]^_<j-A'2 

K[x]l^d K[x]^ ^ K[x]J{a, b, c)^ 0, 

for ly >> 0. As the alternate sum of the dimension of the K-vector spaces is zero and K[x],y/ (a, 5, c)i, 
is for I' >> 0, we have 

= 3(d-j^-l) + i/-/ii-d+l + i/-/i2-d+l + i^ + l 
= d - ^1 - ^2- 

For C{T{x), x" , a;^" ~ I), we have fii + /i2 = 2 n. We are going to show now that in fact /ii = /i2 = n: 
Proposition 2.6. The K[x]-module C{f{x),x'^,x^'^ - 1) has a n-basis. 
Proof. Consider the map 

nAl^i ^ K[a;]3„_i (2.5) 
{p{x),q{x),r{x)) ^ f{x)p{x)+x-^q{x) + (x^" - l)r(x) 

which 3n x Zn matrix is of the form 

/ To -I„ \ 
^ I Ti I„, 1 . (2.6) 



where To,Ti,T2 are the coefficient matrices of {T{x), xT{x), . . . , x"T{x))^ respectively for the list of 
monomials (1, . . . , x""!), (a;", . . . , a;^"-!), (a;^"^ . . . ^ a;^"-!). Notice in particular that T = To + T2 

Reducing the first rows of (To|0| — I„) by the last rows (r2|0|I„), we replace it by the block 
(To + r2|0|0), without changing the rank of S. As T = Tq + T2 is invertible, this shows that the 
matrix S is of rank 3n. Therefore, there is no syzygies in degree n — 1. As the sum 2n = /ii + /i2 and 
A*i < n, H2 < n where ^1, ^2 are the smallest degree of a pair of generators of C{T{x), a;", a;^" — 1) of 
degree < n, we have fJ-i = IJ-2 = n. Thus there exist two linearly independent syzygies (ui, wi, ifi), 
{u2,V2, W2) of degree n, which generate C{T{x), a;", x^" — 1). □ 

A similar result can also be found in [12], but the proof much longer than this one, is based on 
interpolation techniques and explicit computations. Let us now describe how to construct explicitly 
two generators of £(T(x), a;", a:^" — 1) of degree n (see also [T2]). 

As T(a:) is of degree < 2 n — 1 and the map (12. 5p is a surjective function, there exists {u, v, w) G 
K[a;]^_i such that 

f (a;)M(a;) + x"w(x) + (a;^" - 1) w = f (a;)a:", (2.7) 
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we deduce that {ui,vi,wi) — (x" — u, —v, —w) £ C{T{x), x", a;^" — 1). 
As there exists {u',v',w') S K[a::]^_;^ such that 

f{x)u'{x)+x''v'{x) + - = 1 = a;" - (a;^" - 1) (2.8) 

we deduce that (w2, ^2,^2) = (— u', a;" — v', —w' — 1) G £(T'(a;), a;", a;^" - 1). 

Now, the vectors {ui,vi,wi), {u2, V2, W2) of C{T{x), a;", a;^" — 1) are Hnearly independent since by 
construction, the coefficient vectors of a;" in {ui,vi,wi) and (^2,1^2,^2) are respectively (1,0,0) and 
(0,1,0). 

Proposition 2.7. The vectoru is solution of (j2.ip if and only if there existv{x) G K[a:]„_i, w(a;) G 
K[a;]„_i such that 

{u{x),vix),wix)) G £(T'(a;),a;",a;2" - l;g(a;)) 

Proof. The vector u is sohition of (|2.ip if and only if we have 

nE{T{x)u{x)) ^ g{x). 

As u(x) is of degree < n— 1, we deduce from (|2.2p and (|2.3p that there exist polynomial ^(a;) G K[a:]„_2 
and B{x) G K[a:]„_i such that 

r(a;)u(a;) - a;-"+M(a;) - a;"B(a;) = ^(a:). 
By evaluation at the roots lu G il2n, and since w^" = cj" and T{uj) — T(lo) for w G il„, we have 

T(w)u(a;) + uj"v{u;) ~ g(cj), Vlj G il2,i(w), 

with i7(a;) = — a; j4(x) — B{x) of degree < n — 1. We deduce that there exists w{x) G K[a;] such that 

f'(a;)u(a;) + a:"u(a;) + (a;^" - l)u;(a;) = g{x). 

Notice that w{x) is of degree < n — 1, because (x^" — 1) w(a;) is of degree < 3n — 1. 

Conversely, a solution {u{x),v{x),w{x)) G C{T{x), x"^ , x^^ — l;g{x)) nK[a;]^_;^ implies a solution 
{u,v,w) G K"^" of the linear system: 





where S is has the block structure (|2.6p . so that T2 u + w = and Tqu — w = (Tq + T2)u = 5. As we 
have Tq + T2 = T, the vector u is a solution of (|2.ip . which ends the proof of the proposition. □ 



2.2. Euclidean division. As a consequence of proposition [zTol we have the following property: 
Proposition 2.8. Let {{ui,vi,wi),{u2:V2,W2)} an-basis o/£(T(a;), a;", a;^" — 1), the remainder 

/ \ /Ui U2\ 

of the division of a;" g by i Vi V2 is the vector solution given in the proposition (|2.7|) . 
\ g J \wi W2J 

I ° 

Proof. The vector a;"(7 | G jC{T{x), x" , x^^ — l;g) (a particular solution). We divide it by 

\-9. 

'ui U2\ 

Vi V2 we obtain 

.Wi W2I 




(u,w,w) is the remainder of division, thus {u,v,w) G IK[x]^_]^ fl £(r(a;), x", a;^" — l;^)- However 
(ujW, w) is the unique vector € K[x]'^^_i n £(r(x), x", " — l:g) because if there is an other vector 
then their difference is in C{f (x) , x"^ - 1) n which is equal to {(0,0,0)}. □ 

PROBLEM 2.9. Given a matrix and a vector of polynomials ( ^J}^\ ) of degree n, and 



/(^) /'(^) 

Pi^)\ dgg^gg m > n, such that [" ^7 ) invertible; find the remainder of the division of 



^^^y ^ .... ^, 

Pix)\ ( e(a;) e'(x) 



q{x) \f{x) fix) 



Proposition 2.10. The first coordinate of remainder vector of the division of ( _^„^ ) by 



\ , fu u' 



is the polynomial v{x) solution of (12. 

We describe here a generahzed Euchdean division algorithm to solve problem 

Let E{x) = (^gll^l^ °f degree m, B{x) = (^j(^) /'(a^)) '^^S''^*^ n < m. E{x) = B{x)Q{x) + 
R{x) with deg(i?(a;)) < n, and deg((5(a;)) < m — n. Let 2: = ;j 

£;(a;) = B(x)Q(a;) + i?(a;) 

^ Ei-) =B{-)Q{-) + R{-) 

z z z z 

^^^z'^Ei-)^ z"S(i)z™-"g(i) + z"'-"+^z"-^R{-) 
z z z z 

^ E{z) =B{z)Q{z) + z'''-''+^R{z) (2.9) 

with E{z), 13{z),Q(z), R{z) are the polynomials obtained by reversing the order of coefficients of 
polynomials E{z), B{z),Q{z), R{z). 



(1^ ^ = Q{z 



Biz) Biz) 

^Qiz)^EM. modz"^"+i 
Biz) 

— exists because its coefficient of highest degree is invertible. Thus Qiz) is obtained by computing 

Biz) 

the first m — n + 1 coefficients of -J^ . 

Biz) 

To find Wix) = we will use Newton's iteration: Let fiW) = B - W~^. 

Bix) 

f'iWi).iWi+x - Wi) = -Wi-\Wi + 1 - Wi)Wi-^ = fiWi) = B- thus 

Wi+i^2Wi-WiBWi. 

and Wo = B^^ which exists. 

W - Wi+i =W -2Wi + WiBWi 
= Wih - BWi)'' 
= iW - Wi)i3iW - Wi) 

Thus Wiix) = Wix) mod x"^^ for / = 0, . . . , [log(m - n + 1)] . 

Proposition 2.11. We need 0(nlog(n) log(m — n) + mlogm) arithmetic operations to solve 
problem ([^ 

Proof. We must do [log(m — n + 1)] Newton's iteration to obtain the first m ~ n + 1 coeficients 

of — = Wix). And for each iteration we must do 0(77, log n) arithmetic operations (multiplication of 
B 

polynimials of degree n). And then we need 0(771 log m) aritmetic operations to do the multiplication 

E.l.U 

B 
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2.3. Construction of the generators. The canonical basis of K[a;]'^ is denoted by (Ti,o'2,o'3. 
Let pi, p2 the generators of £(T'(a;), x", x^" — 1) of degree n given by 

pi= x'^ai- {u,v,w) = {ui,vi,wi) {21{)') 

P2 = X'^(T2 - [u' ,v' ,w') = {U2,V2W2) 

with {u,v,w), {u',v',w') are the vector given in (|2.7p and (|2.8p . 

We wiU describe here how we compute {ui,vi,wi) and {u2,V2,W2). We wiU give two methods 
to compute them, the second one is the method given in 12J. The first one use the Euchdean gcd 
algorithm: 

We will recal firstly the algebraic and computational properties of the well known extended Eu- 
clidean algorithm (see [T3]): Given p{x),p'{x) two polynomials in degree m and m' respectively, let 

So = 1, si = 0, 

to = 0, ti = 1. 

and define 



Ti+i = Ti^i - qiri, 
Si+i = Sj-i - qiSi, 
ti+i — U-i — Qiti, 

where qi results when the division algorithm is applied to ri_i and r^, i.e. r^-i = g^r^ + r^+i . 
degTi+i < degri for i = 1, . . . , ^ with I is such that r; = 0, therefore ri-i = gcd{p{x) , p' (x)) . 
Proposition 2.12. The following relations hold: 

SiP + tip' = Ti and (si, ti) = 1 for i — 1, . . . , ? 

and 

' degr^+i < degTi, i = l,...,/-l 
degs^+i > degSj and degt^+i > degt^, 
degsj+i = deg(gi.Si) = degw - degr^, 
degti+i = deg{qi.ti) = degu - degr^. 



Proposition 2.13. By applying the Euclidean gcd algorithm in p{x) — x"'~^T andp'{x) — x^""^ 
in degree n — 1 and n ^ 2 we obtain pi and p2 respectively 

Proof. We saw that Tu — g \i and only if there exist A{x) and B{x) such that 

f{x)u{x) + x^'^-^Bix) = x''-^b{x) + A{x) 

with T{x) = a;"~^T(a;) a polynomial of degree < 2n — 2. In (|2.7p and (|2.8p we saw that for g{x) — 1 
((? = ei) and g{x) = x"'T{x) {g = (0, t_„+i, . . . , we obtain a base of £(T'(x), x", a;^" — 1). 

Tui = ei if and only if there exist Ai{x), Bi{x) such that 

T{x)uiix) + x^-^-^Bi{x) ^ + A^{x) (2.11) 

and Tu2 — (0, . . . , i-i)^ if and only if there exist ^2(2^), B2(x) such that 

T{x){u2{x) + x") + x^''-^B2{x) = A2{x) (2.12) 

with degyli(a;) < n — 2 and degyl2(a;) < n — 2. Thus By applying the extended Euclidean algorithm 
in p{x) = x^~^T and p'{x) — x^"~^ until we have degri(a::) = n — 1 and degn+i (x) = n — 2 we obtain 

ui{x) = —si{x), Bi{x) = —ti{x), x"^^ + Ai{x) = —ri{x) 

Ci Ci Ci 
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and 

+ U2{x) = —si+i{x), B2{x) = —ti+i{x), A2[x) = —ri+i{x) 

C2 C2 C2 

with ci and C2 are the highest coefficients of ri(x) and s;+i(x) respectively, in fact: The equation 
(|2.1ip is equivalent to 



n- 1 1 


/ t-n+l 






to 


t-n+1 






■1 










tn-1 


to 






n- 1 1 


\ 


tn~l 


1 


1 



Ul 

Bi 
V / 



1 




since T is invertible then the (2n — 1) x (2n — 1) block at the bottom is invertible and then ui 
and El are unique, therefore wi, Bi and Ai are unique. And, by proposition (|2.12p . degn = n — 1 
{u = ci(a;" + Ai{x)) then degsj+i = (2n— \) — {n~l)=n and degf;+i = (2n — 2) — (n— 1) = n— 1 
thus, by the same proposition, degs; < n — 1 and degi/ <n — 2. Therfore ^s; = ui and = Bi. 
Finaly, Tu = ei if and only if there exist v{x)^ w{x) such that 



T{x) = T^ 



f (a;)u(a;) + + {x^" - l)u'(a;) = 1 

- x2"T- = T + - l)T-thus 

T{x)u{x) + x^vix) + (x^" - l){w{x) + T-{x)u{x)) = 1 



(2.13) 
(2.14) 



ofaother handT{x)u{x)~x-''+^Ai{x) + x''Bi{x) = 1 and 2;-"+Mi(a;) = a;"(a;yli) -a;-"(a;2" - 1)2:^1 
thus 



T{x)u{x) + x"{Bix) - xA{x)) + {x-"" - l)x-"+M(x) = 1 



(2.15) 



By comparing (|2.14p and (|2.15p . and as 1 = a;"x" — (x^" — 1) we have the proposition and we have 



w(a;) 



A{x) — T^{x)u{x) + 1 which is the part of positif degree of —T^{x)u{x) + 1.0 



Remark 2.14. A superfast euclidean gcd algorithm, wich uses no more then 0{nlog^n), is given 
in chapter 11. 

The second methode to compute (ui, vi,wi) and (m2, V2,W2) is given in [12]. We are interested in 
computing the coefficients of ai, (72, the coefficients of 173 correspond to elements in the ideal (x^" — 1) 

a;" - Mo -vo 

u{x) v{x) 
u'{x) v'{x)^ 

A superfast algorithm to compute B{x) is given in [12j . Let us describe how to compute it. 
By evaluation of (|2.10p at the roots ojj G ll2n we deduce that {u{x) v{x))'^ and {u'{x) v'(x))'^ are 
the solution of the following rational interpolation problem: 



and thus can be obtain by reduction of {T{x) x").B{x) by x —1, with B{x) 



f{ujj)u{ujj) + uj'^v{ujj) = 



with 



Un = 1, Vn = 
U'n = 0, < = 1 

Definition 2.15. The r-degree of a vector polynomial 'w{x) = (wiix) W2{x))'^ is defined as 

T — degw{x) :— max{degwi(a:), degW2{x) — t} 
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B{x) is a n— reduced basis of the module of all vector polynomials r{x) G IK[x]^ that satisfy the 
interpolation conditions 

fJriLU,)^0, j=0,...,2n-l 



with fj = 




B{x) is called a r— reduced basis (with t = n) that corresponds to the interpolation data 
i^jJj)^ j = 0, . . . ,2n - 1. 

Definition 2.16. A set of vector polynomial in K[a;]^ is called r-reduced if the r-highest degree 
coefficients are lineary independent. 

Theorem 2.17. Let t = n. Suppose J is a positive integer. Let cti, . . . , cj G K and 0i, . . . , 0j G 
K wich are ^ (0 0)^. Let! <j < J and tj G Z. Suppose that Bj{x) G K[a;]2x2 ^ ^-j -reduced basis 
matrix with basis vectors having tj — degree 6i and 62, respectively, corresponding to the interpolation 
data {((Ti, 4>i); i = 1, . . . 

Let Tj^j := Si — 62. Let Bj^j{x) be a Tj^j-reduced basis matrix corresponding to the interpolation 
data {{ai, Bj{aj)4>i);i = j + I, . . . , J}. 

Then Bj{x) := Bj{x)Bj^j{x) is a tj -reduced basis matrix corresponding to the interpolation data 
{{a^,(l)^);i = 1, . . . , J}. 

Proof. For the proof, see [T2]. □ 

When we apply this theorem for the ujj G il2n as interpolation points, we obtain a superfast 
algorithm {O{n\o^ n)) wich compute 12,^ 
We consider the two following problems: 

3. Bivariate case. Let m G N, m G N. In this section we denote by E' = < « < 

m — 1, < i < n — 1}, and R = K[a:,j/]. We denote by IK[x,?/]m the vector space of bivariate 
polynomials of degree < to in a; and < n in y. 

Notation 3.1. For a block matrix M , of block size n and each block is of size m, we will use the 
following indication : 

M = (Af(ii,i2),(jij2))o<ii ji<m-l = {Mai3)a,l3GE- (3.1) 

(*2, ^2) gives the block's positions, («i, Ji) the position in the blocks. 

PROBLEM 3.2. Given a Toeplitz block Toeplitz matrix T = {ta-i3)aeE.peE G ^"mxmn ^ 
(Taf3)a,P£E With Tap = ta-p) of sizc mn and g = {ga)aGE G K™", find u — iua)aeE such that 

Tu = g (3.2) 



Definition 3.3. We define the following polynomials: 

• T{x,y) := Ujx'y^, 

(j,i)GE-E 
2n-l,2m-l 

• T{x,y) :~ iijx^y-' with 

ti,j si i < m, j < n 

U-2m.,j sii>m,j<n 
1 U.j-2n sii <m,j >n ' 

ti-2m.,j-2n si i > m,i > n 

• u{x,y) -.^ ^ Ui^jx'y' , g{x,y) -.^ ^ gi^x'y^ . 

(ij)eE (i,i)eE 

3.1. Moving hyperplanes. For any non-zero vector of polynomials a = (ai, . . . , a„) G K[x, j/]", 
we denote by C{a.) the set of vectors {hi, . . . , hn) G K[x, j/]" such that 



ti 



J2a^h, = 0. (3.3) 



i=l 



It is a K[a;, yj-module of K[x, y]". 

Proposition 3.4. The vector u is solution of (13. 2p if and only if there exist h2,...,hg £ 
K[x, y]m-i such that y), h2(x, y), . . . , hg[x, y)) belongs to 

n-l 

£(f - ?/",(a;2'" - l)y",y^" - l,x"'{y^" - l),{x^"' - l){y^" - 1)). 

Proof. Let L — {x^^y^^jO < ai < m — 1,0 < a2 < n ~ 1}, and the projection of R on the 
vector space generated by L. By [S], we have 



Tu^ g ^ nE{T{x, y) u{x, y)) = g{x, y) 



(3.4) 



which implies that 

r(x, y)u{x, y) = g{x, y) + y) + x^y-''A2{x, y) + x-^y^A^{x, y) + x'^'y-^A^ix, y) 

+ x^A,{x, y) + x-™A6(x, 2/) + 2/"^7(a;, + y-"^8(a;, y), (3.5) 

where the Ai(x,y) are polynomials of degree at most m — 1 in a; and n — 1 in y. Since = w"'", 
= v~^, T{uj, v) = T{uj, v) for cu G il2m) v & ii2n, we deduce by evaluation at the roots ui G ii2m, 
u e il2 n that 



i?(a;, y) := ^(x, y)u(a;, y) + a;"/i2(a;, y) + y"/i4(a;, y) + x"'-y'^h^{x, y) - g{x, y) G (a; 



2 m 



1,2/'" -1) 



with /i2 = -(^5+^6), hi = -(^7 + ^8), = -{Ai{x,y) + A2{x,y) + A3{x,y) + Ai{x,y)). 

By reduction by the polynomials — 1, y^" — 1, and as R{x,y) is of degree < 3m— 1 in a: and 
< 3n — 1 in y, there exist hs^x, y), he{x, y), . . . , hs{x, y) G K[x, y]m-i such that 

n-l 

fix, y)u{x, y) + x^ h2(x, y) + [x^^ - \)h^[x, y) + y^hi{x, y) + x"'y"h,ix, y) + (3.6) 
(a;2" - l)y-he{x, y) + (y^" - l)/j7(a;, y) + .T™(y2™ - l);i7(x, y) + (.x^" - l)(y2" - l)h^{x, y) = g{x, y). 



Conversely a solution of (|3.6p can be transformed into a solution of (|3.5p , which ends the proof of the 
proposition. □ 

In the following, we are going to denote by T the vector T = {T{x,y), x™, — 1, y", x™ y", (a;^™ — 
l)y",y2" _ l,x'"(y2" - 1), (x^™ - 1) (y2» - l)). 

Proposition 3.5. There is no elements o/K[x,y],„-i in -C(T). 

n-l 

Proof. We consider the map 



^[X,y]m-l K[x,y]3m^l 
n-l 3n-l 

P{x,y) = {pi{x,y), . . . ,pg{x,y)) ^ T.p 



(3.7) 

(3.8) 
(3.9) 



which 9mn x 9mn matrix is of the form 





/ 




E21 


—Ell - 


^Esi 










—Ell 


— E21 


Ell 


— E31 






To 


E2n 


— Ein - 


- Ezn 










— Eln 


— E2n 


Eln 


— Esn 














Ell 


E2I 


—Ell 


+ E31 










5* := 




Ti 








Eln 


E2n 


— Eln 


+ Ezn 






























Ell 


E2I 


-Eii^ 


- E31 




\ 


T2 
















Eln 


E2n - 


-Eln H 


- E^n 



(3.10) 



with Eij is the 3m x mn matrix 6^, (8) Im and is the 3 x n matrix with entries equal zero except 



the («,j)th entrie equal 1. And the matrix 2^1 is the following Qmn x m matrix 



( to 





\ 


tl 


to 





tn-1 




tl to 





tn-l 




t-n+1 











t-n+1 


t-1 




t-n+1 





t-1 


t-n+1 






t-1 



V 



/ 



and i, = 











\ 


ti,i 


tifi 











t,'i t'to 







ti.ni—1 


ti,i 




ti^—7n+l 













ti.—m+ 


1 


ti,-l 




ti,-m+l 







t^,-l 


■ ■ ■ ti—yn + 


1 






ti,-l 




V 







/ 



For the same reasons in the proof of proposition (|2.6p the matrix S is invertible. □ 

Theorem 3.6. For any non-zero vector of polynomials a = (0^)^=1.. G IK[2:,?/]", the ^[x^y\- 

module £(ai, . . . , a„) is free of rank n — \ . 

Proof. Consider first the case where Oi are monomials. 

ai = x""-y^' that are sorted in lexicographic order such that x < y, ai being the biggest and a„ 
the smallest. Then the module of syzygies of a is generated by the S'-polynomials: 

S{ai,aj) = lcm(ai, %■)(— - — ), 



where (CTi)i=i,...^„ is the canonical basis of ]K[a;, y]" [3]. We easily check that S{ai,ak 

lcm(ai,afc) 



S{ai,aj)- 

icm(a"afcj '^("j' "fc) * 7^ J 7^ ^ lcm(ai,aj) divides lcm(ai, afc). Therefore C{a.) is generated by 



lcm(ai,afc) 

lcm(ai ,aj ) 



the S{ai,aj) which are minimal for the division, that is, by 5(0^,0^+1) (for i = 1, . 



1), since 



the monomials are sorted lexicographically. As the syzygies S{ai,ai^i) involve the basis elements 
(Ti,ai+i, they are linearly independent over K[a:, j/], which shows that C{a.) is a free module of rank 
n — 1 and that we have the following resolution: 







[x,y]"-i^K[a;,y]"^(a)^0. 



Suppose now that are general polynomials G K[a;,?/] and let us compute a Grobner basis of a^, 
for a monomial ordering refining the degree [3j. We denote by mi, . . . , the leading terms of the 
polynomials in this Grobner basis, sorted by lexicographic order. 

The previous construction yields a resolution of (mi, . . ., ms): 

^ K[x, y]'-' ^ K[.x, yY ^ (m,)^=i,...,. ^ 0. 
Using [7] (or [3]), this resolution can be deformed into a resolution of (a), of the form 

Q-.K[x,y]P^K[x,yr^{ai)^0, 

which shows that C{a.) is also a free module. Its rankp is necessarily equal to n — 1, since the alternate 
sum of the dimensions of the vector spaces of elements of degree < v in each module of this resolution 
should be 0, for G N. □ 

3.2. Generators and reduction. In this section, we describe an explicit set of generators of 
£(T). The canonical basis of K[2;, y]^ is denoted by ci, . . . , erg. 

First as T'(x, y) is of degree < 2 m — 1 in a; and < 2 n — 1 in y and as the function p.7p in surjective, 
there exists ui, U2 € lK[a;, y]^j_i such that T • ui = T{x, y)x^^, T ■ U2 = T{x, y)y^ ■ Thus, 



pi = x"Vi - ui G C{T), 
P2 =y"fTi G C{T). 
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We also have U3 e K[x, y]m-i, such that T • 1*3 = 1 = a;'" x"' - (a;^™ - 1) = y" - (y^" - l). We 
deduce that 



P4 = y"cr4 - (77 - 1*3 e £(T). 
Finally, we have the obvious relations: 

p5 = y"(T2-(75e/:(T), 

Pe^x"'a4-c75eC{T), 
Ps^y^CTs-Tg + da g£(T). 
Proposition 3.7. The relations pi, . . . ,ps form a basis o/£(T). 

Proof. Let h — {hi, . . . , h^) e 'C(T). By reduction by the previous elements of C{T), we can 
assume that the coefficients hi, /12, /i4, /is are in K[a;, y]m-i- Thus, T{x, y)/ii+x™/i2+y"/i4+x™y"/i5 £ 

ri-l 

(x^ " — 1, y^™ — 1). As this polynomial is of degree < 3 m — 1 in .t and < 3 n — 1 in y, by reduction 
by the polynomials, we deduce that the coefficients /13, h^, . . . ,h<j are in K[a;, y]m-i. By proposition 

n-l 

13.51 there is no non-zero syzygy in K[a;, y]f„_]^. Thus we have h = and every element of £(T) 

Tl-l 

can be reduced to by the previous relations. In other words, pi,. . . ,ps is a generating set of the 
K[x, y]-module C{T). By theorem l3.6| the relations pi cannot be dependent over K[x,y] and thus 
form a basis of >C(T). □ 

3.3. Interpolation. Our aim is now to compute efficiently a system of generators of >C(T). 
More precisely, we are interested in computing the coefficients of cti, CT2, 1T4, tJs of pi,p2, pz. Let us 
call B(x,y) the corresponding coefficient matrix, which is of the form: 



n^Mtn-i (3.11) 

n-l 

Notice that the other coefficients of the relations pi, P2, Pa correspond to elements in the ideal (x^ — 
1, y^" — 1) and thus can be obtained easily by reduction of the entries of {T{x, y), a;™, y", a:™ y^)-B{x, y) 
by the polynomials a;^*" — l,y^" — 1. 

Notice also that the relation p4 can be easily deduced from p^, since we have p3 — x"^a2 + 173 + 
y" (74 — (T7 = p4. Since the other relations pi (for i > A) are explicit and independent of T{x, y), we 
can easily deduce a basis of £(T) from the matrix B{x,y). 

As in C{fT)C\K.[x, y]m-i there is only one element, thus by computing the basis given in proposition 

n-l 

(|3.7p and reducing it we can obtain this element in £(T) n K[x,y]m_i which gives us the solution of 

n-l 

Tu — g. We can give a fast algorithm to do these two step, but a superfast algorithm is not available. 

4. Conclusions. We show in this paper a correlation between the solution of a Toeplitz system 
and the syzygies of polynomials. We generalized this way, and we gave a correlation between the solu- 
tion of a Toeplitz-block- Toeplitz system and the syzygies of bivariate polynomials. In the univariate 
case we could exploit this correlation to give a superfast resolution algorithm. The generalization of 
this technique to the bivariate case is not very clear and it remains an important challenge. 
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